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We derive an expression for the anomalous viscosity in an anisotropically expanding 
quark-gluon-plasma, which arises from interactions of thermal partons with dynamically 
generated color fields. The anomalous viscosity dominates over the collisional viscosity for 
' jf \ \ large velocity gradients or weak coupling. This effect may provide an explanation for the 

apparent "nearly perfect" liquidity of the matter produced in nuclear collisions at RH1C 
without the assumption that it is a strongly coupled state. 

. 81. Introduction 

> . 

o . 

C ■ Measurements of the so-called elliptic flow parameter V2 of hadrons emitted in 

noncentral collisions of heavy nuclei at the Relativistic Heavy Ion Collider (RHIC) 
are in remarkably good agreement with the predictions of ideal relativistic fluid 
dynamics. 1 ) In order to achieve the agreement with the data, the hydrodynamical 
calculations need to assume a rapid initial equilibration within a time 2 ) t\ < 1 fm/c 

i-£h ' and an extremely small ratio of the shear viscosity rj of the fluid to its entropy density 

s. 3 ) In fact, the ratio rj/s cannot be much larger than the conjectured lower bound 
Qh. (47r) _1 for this quantity. 4 ) 

These findings have led to the general conclusion that the matter produced at 
RHIC is strongly interacting with a large opacity, i.e. small mean free path. How- 
ever, this does not necessarily have to be the case: In a previous Letter 5 ) we have 
argued that the unavoidable color field instabilities in a rapidly expanding quark- 
gluon plasma generate an anomalous viscosity, which is smaller than the collisional 
viscosity in the weak coupling limit and for early times after the onset of the longi- 
tudinal expansion. The expression derived in Ref. 5) for the anomalous viscosity, 



VA 1 /(A 2 -1)Tt^ 12 
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(1-1) 



where N c denotes the number of colors and bo is a presently unknown numerical 
coefficient, is parametrically (in the coupling constant g) smaller than the collisional 
viscosity 6 ) 

2° = (1-2) 

where df ~ 5 for three light quark flavors. The anomalous viscosity is generated by 
the diffusive transport of quasithermal partons in the turbulent color fields created 
by Weibel-type instabilities. Such instabilities always arise when the parton momen- 
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turn distribution is anisotropic, 7 )~ 9 ) as is necessarily the case when the quark-gluon 
plasma expands with a preferred direction. The nonlinear interaction of the unstable 
modes has been shown to result in a turbulent cascade of field energy into increasingly 
short wavelength modes, 10 )' 11 ) the characteristic signature of a turbulent plasma. 

The purpose of this article is to present additional details of the calculation 
outlined in our previous Letter. The article is structured as follows. In Section 2 
we discuss the evidence from lattice QCD calculations for a quasi-particle nature of 
the quark-gluon plasma, even in the vicinity of the transition temperature T c . This 
evidence is derived from the off-diagonal flavor structure of the quark susceptibili- 
ties. On the other hand, the RHIC data on jet quenching combined with entropy 
constraints demand a plasma almost entirely composed of gluons. The anomalous 
transport properties owing to the presence of intense coherent color fields can alle- 
viate these constraints. In Section 3 we give a heuristic derivation of the anomalous 
viscosity t/a- In Section 4 we derive the transport equation for quarks and gluons in 
a turbulent plasma. Section 5 contains a brief reminder of the definition of the shear 
viscosity. In Section 6 we present details of the formal derivation of t/a within the 
framework of the theory of transport phenomena in turbulent plasmas. We summa- 
rize our results in Section 7. Two appendices give an alternative derivation of the 
transport equation in the presence of plasma turbulence and provide details on the 
plasma instability structure, respectively. 

§2. The Quark-Gluon Plasma at RHIC 

There is now overwhelming consensus that matter in (approximate) thermal 
equilibrium and with an energy density far in excess of 1 GeV/fm 3 is produced in 
nuclear collisions at RHIC. 1 )' 12 ) The evidence leading to this conclusion rests on two 
key observations. Firstly, all stable hadrons including multistrange baryons are emit- 
ted with chemical equilibrium abundances 13 )' 14 ) and thermal transverse momentum 
spectra that are boosted by a collective transverse ("radial") flow. 15 )' 16 ) Secondly, 
the hadron spectra measured in noncentral collisions show an azimuthal quadrupole 
anisotropy ( "elliptic" flow) characteristic of the fluid dynamical expansion of a fire- 
ball with oval shape. 17 )' 18 ) Plotted as a function of the transverse kinetic energy and 
scaled by the number of valence quarks of the hadron (n = 2 for mesons and n = 3 for 
baryons), the flow parameter v-i exhibits a universal dependence. 19 ) This has been 
interpreted as evidence for a partonic origin of the observed flow pattern, suggesting 
that the transverse expansion of the matter is generated during a phase in which it 
contains independent quasi-particles with the quantum numbers of quarks 20 ) -24 ) . 

The observed magnitude of the elliptic flow requires an early onset of the period 
during which the expansion is governed by fluid dynamics (earlier than 1 fm/c after 
first impact) and nearly ideal fluid properties with a viscosity-to-entropy density 
ratio rj/s <C ]_.25)-28) ^phis resu it ; interpreted in a quasi-particle picture, in turn 
requires that the mean free path in the medium must be extremely short, less than 
the average particle spacing. A short mean free path is also deduced from the strong 
suppression of the emission of hadrons with a transverse momentum pr of several 
GeV/c or more. 29 )' 30 ) This phenomenon, commonly referred to as "jet quenching", 
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is attributed to a large rate of energy loss of the energetic parton in the medium 
before it fragments into the observed hadron. 31 ) -39 ) Since the energy loss is inversely 
proportional to the mean free path of the parton, this observation again requires a 
very short mean free path or, equivalently, a large scattering cross section. 

There are other, purely theoretical reasons for believing that strongly interact- 
ing matter makes a rapid transition in a narrow temperature range around T c 170 
MeV from a hadronic resonance gas to a plasma whose thermodynamic proper- 
ties can be well described by quasiparticle excitations with the quantum numbers 
of quarks and gluons. Phenomenological attempts to describe the thermodynamic 
variables in terms of noninteracting quasiparticles with an effective thermal mass 
work astonishingly well for all temperatures down to T c . 40 )~ 43 ) Rigorous resumma- 
tion techniques using effective quasiparticle propagators for quarks and gluons are 
also quite successful above (1.5 — 2)T C . 44 )' 45 ) 

Even more compelling is the comparison of the temperature dependence of the 
diagonal and off-diagonal quark flavor susceptibilities calculated on the lattice 46 ) 
with expectations from a quasiparticle picture. An especially sensitive quantity is the 
ratio between the baryon number-strangeness correlation (BS) and the strangeness 
fluctuation (.S 2 ), 47 ) which rapidly changes from the behavior characteristic of a 
hadron gas to the temperature independent value of a quasiparticulate quark-gluon 
plasma at T c . 48 ) This result strongly suggests that the quasiparticles that carry 
flavor must have the same quantum numbers as quarks almost immediately above 
T c (for an alternate view, see Ref. 49)). 

However, given the experimental evidence for the low viscosity and strong color 
opacity of the matter, the interpretation of the data in terms of the quasiparti- 
cle picture is problematic: In a quasiparticle picture the conjectured lower bound 
i]/ s > 1/47T corresponds to an extremely short mean free path. Using the relation 
rj ~ nXfp/2 from standard kinetic theory, where n denotes the particle density and 
Xf is the mean free path, as well as the equations p = 3T and s 4n holding for 
massless particle, one finds A/ > (37rT/2) _1 . This implies that at the lower viscosity 
bound the mean free path must not be larger than about half the average distance 
between quasiparticles. This observation is corroborated by simulations of the par- 
ton transport within the framework of the Boltzmann equation with binary elastic 
scattering, which require cross sections up to twenty times larger than expected 
on the basis of perturbative QCD in order to reproduce the elliptic flow data. 50 ) 
Therefore a collisional origin of such a low viscosity is not easily compatible with the 
notion that the medium is composed of well-defined quasiparticles. The quasiparticle 
picture is further cast into doubt by recent calculations of the spectral densities of 
correlators of the stress-energy tensor in an exactly solvable, strongly coupled quan- 
tum field theory (N = 4 super symmetric Yang-Mills theory). 51 )' 52 ) These do not 
reveal peak-like structures that can be attributed to quasiparticle excitations, with 
the sole exception of the hydrodynamical sound mode. 

One possible resolution of the puzzle is to argue 53 ) that the quark-gluon plasma 
near T c is actually strongly coupled. This argument is bolstered by the fact that re- 
sults from solvable strongly coupled gauge theories are in qualitative agreement with 
certain aspects of QCD near T c , such as reduction of the energy density with respect 
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to the ideal gas limit, 54 ) and with experimentally obtained values of transport coef- 
ficients, such as the viscosity, 4 ' 1 ' 55 )' 56 ) energy loss parameter, 57 ) and quark diffusion 
constant. 58 )~ 60 ) It has been argued 49 -*' 61 )' 62 ) that the microscopic structure of such a 
system is dominated by complex bound states of the elementary constituents. Note, 
however, that the diagonal and off-diagonal quark flavor susceptibilities calculated on 
the lattice 46 ) strongly constrain - and in many cases rule out - the existence of bound 
states of the elementary constituents above Tc- In addition, a recent calculation of 
the shear viscosity over entropy ratio in weakly coupled J\f = 4 supersymmetric 
Yang-Mills theory yielded a result many times smaller than the corresponding weak- 
coupling result in QCD. 63 ) This finding therefore may actually suggest that the ratio 
rj/s of QCD near the transition point is several times larger than the lower viscosity 
bound of (47r) _1 for this quantity. 4 ) 

The other possible resolution is that the transport properties of the quark-gluon 
plasma under the experimentally relevant conditions are not governed by collisional 
processes involving perturbative interactions among elementary excitations, but by 
collective phenomena. This situation is not uncommon in plasmas, where coherent 
fields can be spontaneously generated due to instabilities in the field equations in 
the presence of the medium. The occurrence of unstable field modes is a familiar 
phenomenon in electromagnetic plasmas with a charged particle distribution that 
is locally not fully equilibrated. The most relevant of these for our purposes is the 
instability discovered by Weibel, which arises when the momentum distribution of 
charged particles is anisotropic. 64 ) 

It has been known for some time 7 )' 8 )' 65 ) that similar instabilities exist in quark- 
gluon plasmas with a parton momentum distribution that is not in thermal equilib- 
rium. As a result of these instabilities long-range color fields can be excited with 
amplitudes far above the thermal level. The generic nature of such color insta- 
bilities has been recognized only in recent years. 9 )' 66 ) Most work exploring the 
consequences of these instabilities 67 ) -70 ) has been focused on the early stage of the 
collision, when the momentum distribution is highly anisotropic and far from equilib- 
rium. The fields generated by the instabilities drive the parton distribution rapidly 
toward local isotropy and thus toward the hydrodynamical regime. 71 ) However, the 
expansion of the quark-gluon plasma under its own pressure ensures that the matter 
never reaches complete equilibrium, and thus the presence of the color instabilities 
persists even during the period when the matter evolves by viscous hydrodynami- 
cal expansion, although the effect of collisions may weaken the color instabilities. 72 ) 
Since the size of the deviation from kinetic equilibrium is proportional to the viscos- 
ity itself, color instabilities are especially important when the quark-gluon plasma is 
weakly coupled and the collisional shear viscosity is large. By producing an anoma- 
lous contribution to the shear viscosity, extended color fields of a large amplitude 
present a mechanism that may be responsible for the observed small shear viscosity 
of the rapidly expanding quark-gluon plasma studied at RHIC. 5 ) 
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§3. Anomalous Viscosity: Heuristic Derivation 



Anomalous contributions to transport coefficients caused by the action of tur- 
bulent fields are well known in electromagnetic plasmas. As we already mentioned, 
the condition for the spontaneous formation of extended electromagnetic fields in a 
plasma is the existence of instabilities in the field equations due to the interaction 
with the charged particles. This condition is satisfied in electromagnetic plasmas 
with an anisotropic momentum distribution. 64 ) The treatment of the consequences 
for transport processes is based on the formalism of particle propagation in turbulent 
plasmas originally developed by Dupree. 73 )' 74 ) The term plasma turbulence refers 
to the spectral distribution of the field excitations, which follows a power law, in 
analogy to the Kolmogorov spectrum of vortex excitations in a fluid with fully de- 
veloped turbulence. Such plasmas are characterized by strongly excited random field 
modes in certain regimes of instability, which coherently deflect the charged parti- 
cles and thus reduce the effective mean free path. The scattering by turbulent fields 
in electromagnetic plasmas is known to greatly reduce the heat conductivity 75 )' 76 ) 
and shear viscosity 77 )' 78 ) of the plasma and to increase the energy loss of charged 
particles propagating through it. 79 ) 

Following Abe and Niu 78 ) the contribution from turbulent fields to plasma trans- 
port coefficients is called anomalous. We shall see below that this designation is 
justified by the fact that the viscous corrections to the hydrodynamic energy-stress 
tensor due to the turbulent fields are nonlinear in the velocity gradient, but exhibit a 
sublinear dependence. As a consequence, the viscous effects in a quark-gluon plasma 
with a large imprinted velocity gradient are much smaller than expected from the 
usual linear theory. As the formation mechanism of quasi-thermal QCD matter in 
relativistic heavy ion collisions naturally leads to a large longitudinal velocity gradi- 
ent in the direction of the beam axis, the relevance of such anomalous contributions 
to various transport coefficients, including the shear viscosity, is not unexpected, 
especially during the early stages of the expansion. 

Before describing the derivation of the anomalous shear viscosity ??a of an ex- 
panding quark-gluon plasma in detail, it is useful to give a heuristic argument for 
its dependence on the amplitude of turbulent plasma fields. This argument will also 
elucidate the reason for the dominance of the anomalous viscosity at weak coupling. 
According to classical transport theory, the shear viscosity is given by an expression 
of the form 80 ) 



where n denotes the particle density, p is the thermal momentum, and Af the mean 
free path. For a weakly coupled quark-gluon plasma, n « 5T 3 and p f» 3T. The 
mean free path depends on the mechanism under consideration. The collisional shear 
viscosity rjc is obtained by expressing the mean free path in terms of the transport 
cross section 



n 



(3-1) 



\f> = (na tl )- 1 . 
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Using the perturbative QCD expression 

5g 4 y/A^ 
Ap 

for the transport cross section in a quark-gluon plasma yields the result 



T 187TS 
77(7 ~ ~ 7= , (3-4) 

<rtr 25g 4 \n(V^/g) V ; 

where we used the relation s ~ An valid for ulrarelativistic particles. This agrees 
parametrically with the result (1-2) for the collisional shear viscosity in leading log- 
arithmic approximation. 6 )*^ 

The anomalous viscosity is determined by the same relation (3T) for 77, but the 
mean free path is now obtained by counting the number of color field domains a 
thermal parton has to traverse in order to "forget" its original direction of motion. 
If we denote the field strength generically by B a (a denotes the color index), a 
single coherent domain of size r m causes a momentum deflection of the order of 
Ap ~ gQ a B a r m , where Q a is the color charge of the parton. If different field domains 
are uncorrelated, the mean free path due to the action of the turbulent fields is given 

Ai*>= M (p/^-^_. (3-5) 

The anomalous shear viscosity thus takes the form: 

np 3 9sT 3 
VA ~ 3g 2 Q 2 (B 2 )r m ~ 4,? 2 Q 2 </3 2 >r m ' W 

which agrees with expression (16) of Ref. 5), if we identify r m with the memory time 
r m for relativistic partons. 

The argument now comes down to an estimate of the average field intensity (B 2 ) 
and size r m of a domain. We first note that the size is given by the characteristic 
wave length of the unstable field modes. Near thermal equilibrium, the parameter 
describing the influence of hard thermal partons on the soft color field modes is the 
color-electric screening mass ttjd ~ gT. Introducing a dimensionless parameter £ 
for the magnitude of the momentum space anisotropy, 9 ) the wave vector domain 
of unstable modes is k 2 < ^m 2 ^ (see Appendix B). Thus r m ~ ^^(gT)" 1 . The 
exponential growth of the unstable soft field modes is saturated, when the nonlin- 
earities in the Yang-Mills equation become of the same order as the gradient term: 
g\A\ ~ k, which implies that the field energy in the unstable mode is of the order of 
g 2 (B 2 } ~ k A ~ £ 2?Ti D- The denominator in (3-6) thus has the characteristic size, at 
saturation: 

ff 2 Q 2 (/3 2 )r m ~e 3/2 ^~e 3/2 (5T) 3 . (3-7) 



The logarithm in (3-3) is the weak coupling limit of the function 80 - 1 I(a s ) = (2a s + 1) ln(l + 
aj 1 ) — 2, which never becomes negative. In fact, a 2 s I{ct s ) — > 1/6 for a s ^> 1, which suggests that 
rj/s approaches a constant in the strong coupling limit. Of course, the derivation of (3-3) becomes 
invalid in this limit. 
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Inserting this result into the expression (3-6) gives the following relation for the 
anomalous viscosity: 

We conclude that the anomalous viscosity will be smaller than the collisional viscos- 
ity, if the coupling constant g is sufficiently small and the anisotropy parameter £ is 
sufficiently large. 

§4. Diffusive Vlasov-Boltzmann Equation 

In order to present a more rigorous derivation of the anomalous viscosity of 
a turbulent quark-gluon plasma, we first need to derive the appropriate transport 
equation. The propagation of quasi-thermal partons in the presence of soft, locally 
coherent color fields and hard collisions among the partons is described by a Vlasov- 
Boltzmann equation: 81 ^ 

v^f(r, P, t) + gF a ■ V p f a (r, p, t) + C[f] = . (4-1) 

Here /(r,p,t) denotes the usual parton distribution in phase space, which sums 
over all parton colors. f a (r, p, t) denotes the color octet distribution function, which 
weights each parton with its color charge Q a . Both, / and f a can be defined in the 
semiclassical formalism 82 ) as the moments of the distribution function f(r,p,Q,t) 
in an extended phase space that includes the color sector: 

/(r,p,t) = J dQf(r,p,Q,t), (4-2a) 

f a (r, p,t) = J dQ Q a f(r, p, Q, t) . (4-2b) 

In (4-1) we have used the covariant notation = p^/p° with p^ = (E p , p). v = p/E p 
denotes the velocity of a parton with momentum p and energy E p . Furthermore, 

F a = £ a + vxB a (4-3) 

denotes the color Lorentz force, and C[f] stands for the collision term. (We will 
specify C[f] later.) 

The color octet distribution f a satisfies a transport equation of its own, which 
couples it to phase space distributions of even higher color-SU(3) representations. 
In the vicinity of the equilibrium distribution, however, it makes sense to truncate 
the hierarchy at the level of the color singlet and octet distributions and to retain 
only the lowest terms in the gradient expansion. We also note that the color octet 
distribution function vanishes in equilibrium, f§ = 0, implying that f a is at least of 
first order in the perturbation. 

The transport equation for f a then reads: 81 )' 82 ) 

+ 9fa bc A b ^f + j^jF a • V p f + C a [f, fl = 0, (4-4) 



c 
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where C2 denotes the quadratic Casimir invariant of the color representation of the 
thermal partons. 

Before we linearize the Vlasov-Boltzmann equation (4-1), we must rewrite it in 
a form applicable to the case of a turbulent quark-gluon plasma. In order to do so, 
we need to make additional assumptions about the field distribution in the Vlasov 
force term. Based on the arguments presented in Section II, we shall assume that the 
color field is turbulent, i. e. random with a certain spatial and temporal correlation 
structure for fields at different space-time points. This assumption will allow us 
to rewrite the force term involving the color octet distribution function f a into a 
dissipative (Langevin) term acting on the color singlet distribution /. 

We shall derive the diffusion term in two different ways. In this Section, we will 
use linear response theory to calculate f a from the color singlet distribution for given 
color fields. The diffusion term is then obtained after substituting the result into 
the force term in (4-1) and choosing appropriate correlation functions for the color 
fields. In Appendix A, we shall present another derivation based on the extended 
distribution function f(r,p,Q,t). There we also discuss the difference between our 
approach and the one due to Dupree 73 ) for turbulent abelian plasmas. 

In order to resolve eq. (4-4) for the color octet distribution, we Fourier transform 
the dependence on the space-time coordinate x M = (t, r): 

/d A k 
—e- ik -*r(p,k). (4-5) 

We allow for an arbitrary particle distribution in momentum space, but neglect any 
space-time dependence of the singlet distribution /(p). Ignoring the collision term 
and, for the moment, the gauge connection associated with the space-time derivative, 
the solution of eq. (4-4) is given by: 83 ) 

(p, fc) = -ig^-^v ■ k + ie^F^k) ■ V p /(p) , (4-6) 

where v ■ k = v^k^ = k° — v • k. The gauge connection has the effect of adding a 
path-ordered factor 

U ac (x, x') = P exp (- £ f abc Al dx^j , (4-7) 

which parallel transports the gauge fields from x' to x. Returning to coordinate 
space, (4-6) then takes the form: 

/"(P.«) = -^ /(0/ d A X 'U ab {x,x')^^F\x').V p f( P ). (4-8) 

Inserting this solution for f a into eq. (4T) for the singlet distribution function /, the 
Vlasov force term takes the following form: 

g F a (x)-V p r(p,x) = -^ l F a (x)-V p J ^0 J d*x'U ab (x,x') 
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We now need to invoke the argument that the soft color fields are turbulent and 
that their action on the quasi-thermal partons in (4-9) can be described by taking 
an ensemble average, which can be factorized in the form 

{F-(x)U ab (x,x')F^x')f(p)) = <i?(z)LUz,z')^V)>/(p). (4-10) 

where / = (/). We furthermore assume that the correlation functions of fields at 
different space-time points x and x' depend only on \x — x'\ and fall off rapidly with 
correlation time r and correlation length a: 

(£^x)U ab (x,x')£^x')) = (\t-t'\) ^ (|x-x'|) , (4-lla) 

{B*(x)U ab (x,x')B){x>)) = (B?B«)<P^ (\t-t'\) &r g) (|x-x'|) . (4-llb) 
Examples satisfying these assumptions are the Gaussian correlators 

# (ei/ mag ) (| t _ f|) = exp [(* _ t') 2 /2r c 2 1/mag ] , (442a) 

^el/mag) (| x _ x '|) = exp [ (x _ X ') 2 /2a c 2 1/mag ] . (4-12b) 

Finally, we assume that the color-electric and -magnetic fields are uncorrelated: 
(£?(x)U ab (x,x')B b j (x'))=0. 

The reality of the correlation functions (4-12) and their symmetry with respect to 
exchange of the two space-time arguments express the chaotic nature of the plasma. 
It is here where the reversibility of the mean field dynamics is explicitly broken 
and the dissipative nature of the turbulent plasma is introduced. The existence of 
plasma instabilities and the associated exponential growth of the unstable modes, 
which correspond to the presence of positive Lyapunov exponents in the coupled 
field-particle system, forms the essential physical basis for our argument. 

Since the right-hand side of (4-9) must be real, under the conditions outlined 
above only the imaginary part of the propagator 

Im = -ir6(vk) (4-13) 

v ■ k + ie 

can contribute. Performing first the integral over k°, then the integral over k and 
finally the integral over x', we obtain: 



/(p) 



= -V p -I>(p)-V p /(p). (4-14) 

Here we made use of the identity (v x B) • V p = —B • (v x V p ) and introduced the 
memory time 



el/mag 

2 



1 r°° 

- / dt'$f!™3 (\t - t'\) c^/mag) (| v ( t _ . ( 4 .15) 

2 J-oo 



The precise matrix structure of the correlators (4-11) may differ between the 
various stages of a relativistic heavy ion collision. Initially, color-magnetic field 
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modes exhibit the strongest growth rate, 84 ) but in the later turbulent steady-state 
regime all spatial components of the color field modes populated by the instabilities 
are expected to be of comparable strength. 10 )> n ) For completeness, we shall discuss 
both scenarios. 

When we consider only the color-magnetic fields initially generated by the plasma 
instability, which point in a transverse direction, we can expect the ensemble average 
to be transverse in space: 

(A) {B?B°) = \{5 t3 - 6 iZ 5 JZ )(B 2 ) ; {S?S?) « . (4-16) 

We shall call this case Scenario A. In the late turbulent phase all spatial components 
of the field correlators are of approximately equal size, and we write: 

(B) <B?B?) = 1 -5 %3 (B 2 ) ■ = ^(£ 2 ) . (4-17) 

We shall call this case Scenario B. Employing the notation — ip x V p = L/ p ) for the 
generator of rotations in momentum space, we can write the diffusive term in the 
transverse color-magnetic field dominated Scenario A as 

V P • D • V p = - 2{N f^\ )E 2 ( g2 ) C ag [(L (p) ) 2 - (4 P) ) 2 ] , (448) 

and in the Scenario B with isotropically turbulent color-electric and color-magnetic 
field excitations as 



9 2 C 2 



(£ 2 )T*Vl-(B 2 )T^ {LiP))2 



V \~ / 'm F 2 



(4-19) 



Note that the operator V 2 associated with random color-electric fields yields a 
nonvanishing contribution when acting on the equilibrated momentum distribution 
/o(p)- This is different for the operator associated with random color- magnetic 
fields, (L/ p )) 2 , which only yields a nonvanishing contribution when acting on the 
anisotropic part of the momentum distribution given by /i(p). This is not surpris- 
ing, because randomly distributed electric fields are well known to lead to an increase 
in the average energy of the plasma particles, 94 ) corresponding to a heating of the 
plasma. In contrast, color-magnetic fields only contribute to the isotropization of 
the momentum distribution and do not cause plasma heating. 

§5. Linear Response Theory 

We now have motivated the replacement of the Vlasov-Boltzmann equation (4T) 
by Dupree's ensemble averaged, diffusive Vlasov-Boltzmann equation 

„ d 



f-V p -D.V p f + {C[f}) = 0. (5-1) 



We next need to expand this equation up to linear terms in the ensemble averaged 
perturbation Sf of the parton distribution. The drift term gives rise to gradients 
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of the collective variables T and which are considered to be of the same size as 
terms linear in Sf; the collision term vanishes at equilibrium and contributes only 
at first order in Sf; and the Vlasov term will require special consideration. 

The general linear response (Chapman-Enskog) formalism assumes a small per- 
turbation of the thermal equilibrium distribution 

/o(p) = (e^Tl) _1 , (5-2) 

where (3 = 1/T denotes the inverse temperature and the upper/lower sign applies to 
bosons/fermions, respectively. Here denotes the four-velocity of the equilibrated 
medium and u • p = u^p^. Using the relation T dfo/d(u • p) = —fo(l ± / ), one 
expresses the perturbed distribution function in the form: 

/(P, r) = /o(p) + Sf(p, r) = / (p)[l + h(p, r)(l ± / (p))], (5-3) 

where the function /i(p, r) can be considered as (minus) the change in the argument 
of the equilibrium distribution: 85 ) 

/(p,r)«(e^ r 'Tl)" 1 . (5-4) 

While the relativistic dissipative hydrodynamics can be cast in a manifestly 
covariant formalism, 86 )' 87 ) this is not necessary for the derivation of the transport 
coefficients, which are Lorentz invariants. It is thus convenient 88 ^ to work in the 
local rest frame of the fluid, i. e. in the frame where u(x) = 0,u° = 1 and thus 
u ■ p = p° = E p at the considered space-time point x. This choice also implies the 
relation dvP/dx 11 = 0. 

In the derivation of transport coefficients one assumes that the particle distri- 
bution is slowly varying, so that the deviation from the equilibrium distribution is 
homogeneous in space and proportional to gradients of the equilibrium parameters. 
For the shear viscosity one uses the local perturbation*) 

/i(p,r) = -^r^PiPjiVujij, (5-5) 

where A(p) is a scalar function of the momentum p, which measures the magnitude 
of the deviation from equilibrium, and (Vu)^- denotes the traceless symmetrized 
velocity gradient: 

(Vu)ij = -(ViUj + VjUi) - -5ijV ■ u, (5-6) 
For Bjorken's 89 ) boost invariant flow field u z = z/t we have 

(Vu)^- = -Uiag(-l,-l,2). (5-7) 
or 



*-* Our definition makes Zi(p) dimensionless. To compare with Ref. 6), identify A(p) = 
— x(p)E P /p 2 . In order to compare with Ref. 85), identify = x(p)/T and A(p) = —g(p)T/E v . 

Also note that (Vu)ij projects out the traceless part of PiPj. 
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The connection to the shear viscosity is made by comparing the microscopic 
definition of the stress tensor 

T -* = / ' iSk; Kn,{v - r> (5 ' 8) 

with the macroscopic definition of the viscous stress: 

T ik = + 5T tk = P5 tk + eu t u k - 2r){Vu) ik - (5 ik V • u, (5-9) 

where e and P are the equilibrium energy density and pressure, and where r\ and £ 
denote the shear and bulk viscosity, respectively. 

Combining equations (5-3), (5-5) and (5-8) one finds 

(2tt) 3 £ P&Mp)fo(p)0- ± /o(p)) 

= -r](Vu) mn [6 iTn Skn + S in 6 km + S ik S mn ], (5-10) 
with the shear viscosity coefficient 

d 3 p p 



T] ~ IhTj (2vr)3 El A[P) dE p {b U) 



For the special case (5-7) of boost-invariant longitudinal flow, the perturbation 
of the equilibrium distribution takes the form 

A<p)~f£|V«| (»«-£), (5-12) 
with (note the normalization): 

1/2 1 



^(Vu)ij(Vu)ji 



T 



(5-13) 



and the anisotropy of the stress tensor (5-10) is 3 )' 90 ) 



2ST XX = 25T yy = -5T ZZ = ^ r?|V«| . (5-14) 
§6. Shear Viscosity 



6.1. Drift Term 

We begin with the evaluation of the drift term. The dominant contribution 
comes from the r-dependence of the local equilibrium distribution. Because the 
form of this term does not depend on ensemble averaging, we shall omit the "bar" 
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symbol in this section. Using a dot to denote the partial time derivative, one finds 
in the local rest frame (recall that v = p/E p ): 

„ d 



/o(p) = -/o(l ± /o) (3E p - (3u ■ p + v • V/3 - 0v- V(u • p) . (6-1) 



This expression can be simplified with the help of the energy-momentum conservation 
law in the presence of color fields: 

-g^r = F aiXV H • (6-2) 

With the help of the expression for T^ v in terms of the momentum space distribution 
function, we find: 

dT» v _ f d 3 p M v dfo 
dx» J (2tt) 3 E p P P dx» 

= dx~^ -a^'^J pfo[p) - (6 ' 3) 

Using the expression (4-8) for the color octet distribution induced by the color field 
and following the same arguments that led to the relation (4-14), we obtain for the 
ensemble average of the right-hand side of (6-2): 



(F a ^(x)f u (x))=9 j J^e^^^^ 



g 2 C 2 f d 3 p 



iVf - 1 J (2tt) 3 E 1 

The time-like component {y = 0) of the expressions (6-3) and (6-4) are easily evalu- 
ated to yield the relation 

4- (v ' u < = <£ "' j ° > = f < ^-- (6 ' 5) 

where mo is the Debye screening mass defined as*-* 

, _ i-C, [ tPp „ „ ,-, ,_/JV c , ATA w 



"* " "SRI J (^)% p ■ V ' /(p) " If + -f) ^' • < M > 

and the final expression holds for a noninteracting plasma of massless quarks and 
gluons in thermal equilibrium. The space-like components yield, after some calcula- 
tions: 

(V/?-/M)^ = -(£ a f° + j a x B a ) = - g J ^(F a r(p)) = 0. (6-7) 



In the notation of Ref. 5) the coefficient of <? 2 T 2 is given by the expression ^$rr^2 
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The right-hand side vanishes by virtue of (4-8), because the momentum integral 
reduces to a surface term at infinity. 

We can now use the relations (6-5) and (6-7) to eliminate the time derivatives (5 
and u from the drift term (6-1), obtaining: 



: ' / (p) = /o(l±/o) 



(p-V)(u-p) ml(£ 2 )r^E p 
E P T 3T 2 (de/dT) 
dP/dT 



T(de/dT) 



E P (V ■ u) 



(6-8) 



The first term in the square brackets can be expressed in terms of the traceless 
velocity gradient (5-6) and a term which can be combined with the last term in the 
brackets, yielding 

For a gas of massless noninteracting partons, e = 3P and p 2 = E p , causing the term 
proportional to the divergence of the collective velocity to vanish, in accordance with 
the expectation that the bulk viscosity of a scale invariant system must be zero. Thus 
our final expression for the drift term is: 



^/o(P) = /o(l ± /o) 



E p T [ h 3T2(de/dT) 



(6-10) 



Equation (6-5) describes the heating of the plasma by the coherent color field. 
The coefficient of (£ 2 ), 

°A = \mlTt, (6-11) 

is the effective color conductivity of the turbulent plasma (see Appendix B of Ref. 83)). 
This expression for the conductivity differs from the collisional conductivity obtained 
for a plasma in equilibrium in the absence of coherent color fields. 91 )> 92 ) In anal- 
ogy to the anomalous shear viscosity t/a of the turbulent plasma, which is the main 
object of this article, <ta may be called the anomalous color conductivity. 

6.2. Force Term 

Next we come to the diffusive Vlasov force term in (5-1). Because they are 
associated with different operators in momentum space (see eq. (4-19)), the color- 
electric and color-magnetic contributions to the diffusion term need to be considered 
separately. Since /o depends solely on E p = |p| and hence L^/q = 0, the color- 
magnetic contribution to the diffusion term affects only the deviation of the particle 
distribution from equilibrium. On the other hand, the color-electric contribution to 
the diffusion term affects the full momentum space distribution, because V 2 fo / 0. 
The difference is easy to understand: color- magnetic fields only change the direction 
of the momentum of a parton, leading to a rearrangement of particles within the 
equilibrium distribution, but not to a modification of the distribution itself. Color- 
electric fields, on the other hand, accelerate partons and thus lead, on average, to a 
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heating of the thermal distribution. In the macroscopic formulation, this difference 
is related to the fact that color-electric fields induce a color current in the parton dis- 
tribution, which interacts dissipatively with the color-electric field. We have already 
discussed this effect at the end of the previous section. 

We first consider the diffusion operator (4-18) for Scenario A. Since the angu- 
lar dependence of the perturbation /i(p) in (5-5) has the form of a quadrupole in 
momentum space, it is an eigenfunction of the magnetic diffusion operator. For the 
perturbation (5-12) associated with the Bjorken flow, /i(p) ~ Y2o(p)> implying that 



(Z>)) 2 - (L^) 2 1 /x(p) = \(L^) 2 ] f^p) = 6/i(p). 



(6-12) 



The color-magnetic part of the diffusion term in Scenario A thus takes the form 



V p • ■ V p /(p) = 



3C 2 A g 2 (B 2 }r l 
(N 2 



mag 



(6-13) 



The diffusion term for scenario B yields (with p 

C 2 g 2 (S 2 )r^ d 2 



IPl): 



V p -D- V p /(p) 



[Pf(p)] 



S(N 2 -l)p dp 2 

2C 2 Ag 2 f(£ 2 



2\ mag 
' m 



E 2 

p 



[N 2 - l)E p T 2 V P 

x/o(l ± fo)PiPj{Vu)ij 



(6-14) 



The rotationally symmetric part of the first term describes the heating of the parton 
distribution by the turbulent electric fields discussed in conjunction with the anoma- 
lous color conductivity (6-11). The anisotropic part of the first term, together with 
the second term describes the angular diffusion of the parton distribution, which 
leads to an anomalous viscosity. These terms have the same structure as the force 
term in Scenario A except that color-electric as well as color-magnetic fields con- 
tribute. 

6.3. Collision Term 

In the evaluation of the collision term 



C[f]( P ) 



1 



d 3 kd 3 p'd 3 k' 



T 8\p + k - p' - k') 



IE P J (2vfSE k E> p E' k 

x(/(p)/(k)[l±/(p')][l±/(k')] 

-/(pO/(k')[l±/(p)][l±/(k)]) 



J2M a (p,k,p',k': 



(6-15) 



we follow Arnold et al9 Note that (6-15) involves an implicit summation over 
parton flavors and helicities. Since the collision term vanishes at equilibrium owing 
to detailed balance, the leading contribution is linear in /i(p). In first approximation, 
the collision term thus gives rise to a linear integral operator of the form 

/o9k 
^Axi2 Wo( P )/o(k) [/i(p) + /i(k) - /!(p') - ji(k')] , (6-16) 
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where 



d(Jl2 = 8^li P ^E'f iP + k ~ p/ ~ k,) 



[l±/ (p')][l±/o(k')] (6-17) 



denotes the differential cross section for the scattering process p,k — > p',k' and 
v re \ is the relative velocity. We also note that ensemble average of the collision 
term required in the diffusive Vlasov-Boltzmann equation (5-1) simply translates 
into the averaged distribution function ]\ in the linearized collision term (6-16): 
(I[fi]) = I[fi]- 

For the leading logarithmic limit of the collisional viscosity, it is sufficient to 
use the contributions to the scattering matrix element M a with the highest infrared 
divergence. Using Mandelstam variables s, t, u, the squared matrix element for one- 
gluon exchange processes (averaged over initial-state and summed over final-state 
quantum numbers) is 

\Mw\ 2 = j^£cP($\ (6-18) 

where a, b denote the quantum numbers of the scattering partons, and are the 

quadratic Casimir operators for their color multiplet. The squared matrix element 
for the quark annihilation process is 

Finally, the Compton scattering process doubles the contribution from quark anni- 
hilation. 6 ) 

6.4. Anomalous Viscosity 

We are now ready to calculate the coefficient of shear viscosity. We begin by 
ignoring the collision term and calculate only the contribution of the diffusive Vlasov 
term, i. e. the anomalous shear viscosity. Equating the first term in (6-10) with the 
right-hand side of (6-13) for Scenario A we obtain 

(N? - \)ElT 

We note that the diffusive Vlasov equations for quarks and gluons decouple in the 
absence of collisions, causing the function A(p) to take different values for quarks 
and gluons. Inserting (6-20) into the relation (5-11) between rj and A and assuming 
massless partons, the desired expression for the anomalous shear viscosity due to the 
action of the coherent color-magnetic fields on massless partons is found to be: 



poo 

VA = , r _ ^%,^ mag / dppVo(p) , (6-21) 

Jo 



157r 2 C 2 <7 2 (fi 2 )rr 



where a sum over parton species and helicities is implied. After performing the 
momentum space integration, one obtains the following anomalous viscosities for 
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gluons and quarks: 

{9) _ igqgKjvg - 1) 2 t 6 

{q) _ 62C(6)iV c 2 iV / T« 

These results differ by a numerical factor from the one obtained previously 5 ) by 
approximating A(p) as a constant. We emphasize again that, given an ensemble of 
color fields, gluons and each flavor of quarks generate their separate contribution to 
the shear viscosity. 

The results obtained for Scenario B have a similar, but somewhat more compli- 
cated form, because the algebraic equation for A{p) is replaced with a second-order 
differential equation. We will not discuss this case further here and come back to it 
in the context of the calculation of the complete shear viscosity. 

However, we note that the field ensemble is sensitive to the total value of tja, 
because the field instabilities are driven by the overall anisotropy of the parton 
distribution. Possessing a smaller color charge than gluons, quarks develop a larger 
momentum space anisotropy in an expanding quark-gluon plasma, neglecting the 
effect of collisions. On the other hand, quarks contribute with a smaller weight than 
gluons to the color polarization tensor, which drives the instability of soft modes of 
the color field. Thus, although A(p) and tja are inversely proportional to C2, the 
contribution of gluons and quarks to the color field instabilities, which is proportional 
to C2?7a, is independent of the magnitude of their color charge. 

6.5. Complete Shear Viscosity 

The full linearized diffusive Vlasov-Boltzmann equation (5-1) constitutes a lin- 
ear integral equation for the scalar function A(p) characterizing the deviation of 
the momentum distribution from equilibrium. An exact solution of this equation 
requires numerical methods. The standard approach 85 )' 93 ) makes use of the fact 
that the kernel of the linearized collision operator is self-adjoint with respect to an 
appropriately chosen scalar product and has non- negative eigenvalues. As a conse- 
quence, the solution of the Boltzmann equation coincides with the minimum of a 
quadratic functional W[/i] and can thus be obtained from a variational principle. 
The variational principle can be cast into the form: 



^ME> + i (_ Vp . D ■ V p 5f(p) + JIAKp)) 



= mm. 



(6-23) 

The minimum of (6-23) defines the optimal solution /i(p) or A(p) of the transport 
equation (5-1). 

The optimal function A(p) can be determined by means of the variational 
method after an expansion in a complete set of orthogonal functions. A good (up 
to a few percent) approximation can be obtained by choosing the one-parameter 
function A(p) = A|p|/T, where A is a constant. Taking the appropriate moment 
of the linearized transport equation results in an algebraic equation for A, which 
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can be solved to obtain an approximate analytic expression for the shear viscosity. 
Because the mean free paths of quarks and gluons are different, we need to introduce 
different parameters for quarks (A q ) and gluons (A g ). 

We now evaluate the momentum integrals in (6-23) for massless quarks and 
gluons. For the drift term we obtain: 



M-'dI/i 



1 f d 3 p 



Tps j j^^A(p)[p iPj (Vu) tJ Yf (p)(l±f (p))v.V r f 



(p) 



2\Vu\ 2 
' 9tt 2 T 3 
32jVuf 
3vr 2 



roo 

/ dpp*Af (p) 
Jo 



C(5)T 2 



15 



{N£ - l)Ag + —N c NfA q 



(6-24) 



The diffusive Vlasov term yields (in Scenario A): 



(A) 



V 



[h\ = -\j 



^3/ 1 (p)V p -L>-V p <5/(p) 



3C 2 g 2 (B 2 )r^ f d*p 2 , un2fnn + fn , 

2T 4( jV 2 _ fj J ^ WiPjiVuhj} /o(p)(l±/o(p)) 



15T 5 iV 2 - 1 



2/r2\ mag ,>oo 



16|Vu| ; 
5vr 2 T 



2 /R 2 W ma g 



dpp 3 A 2 f (p) 
N c A 2 g + 



(6-25) 



For Scenario B, the diffusive Vlasov term receives contributions from color-electric 
as well as color-magnetic fields: 



i = 



d 3 p 
(2^)3 



/i(p)V p -Z>-V p <S/(p) 

d 3 p /<£ 2 )r^ , (B 2 )t, 



C 2 g 2 f d*p f(£ 2 ) 
T\N 2 -I) J (2tt) 3 E 2 V p 



+ 



mag 
m 



E 2 

p 



xA(p) 2 [piPjiVu)^} 2 / (p)(l ± / (p)) 



<i 3 p 



6T 4 (iV 2 - 1) J {2irfp 3 E p 



-A{p) IpiPjiVujijf 



8\Vu\ 2 C 2 g 
45T 5 (iV 2 - 



/ 

dp 2 

2 poo 

T) (2(£ 2 ) rt + (B 2 ) rr s ) J dpp 3 A 2 f (p) 



p 3 A(p) 

E n 



/o(p)(1±/o(p)) 



= ^5^ c(%2(2(£:2)r - + ^ 2)rrs) 



iV c A 2 + -iV>A 2 



. (6-26) 
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This result differs from the one obtained for Scenario A only by the substitution 

2 



mag 
m 



(2(S 2 )r£ + (B 2 }rr e ) ■ 



(6-27) 



Finally, we simply state the result for the collision term in the leading logarithmic 
approximation and refer to Arnold et a/. 6 ) for details of the calculation: 



W c [h 



d 3 p 
(2^)3 

ir\Vu\ 2 T 2 
45 



/i(pWi](p) 



(A^-VlncT 1 

7T 2 N 



l(2N c + N f ) ^N c A 2 g + 7 -N f A^ 



+ 



256N, 



Un 2 - l)(A g - A q f 



(6-28) 



Note that our expression differs from that of Ref. 6) by an overall factor 4| X7u\ 2 / (45T) 
owing to the different definition of the expectation value (6-23). 

In order to calculate the shear viscosity, we have to minimize (6-23) with respect 
to A g and A q , and then insert the obtained values into the expression 



V 



24C(5)T 3 
3tt 2 



(N 2 - l)A g + ^-N c N f A q 



(6-29) 



obtained by performing the momentum integral in the formula (5T1) for the shear 
viscosity. The mininization results in two linear equations for A g and A q . We state 
the expression for Scenario A: 



32C(5) 
3vr 2 



7V 2 -1 
fN c N f 



32C(4) g 2 {B 2 )r i 



mag 
m 



5vr 2 
^(iV 2 



45 



T 3 



1 4 1 -1 
-9 In g 



n 2 N f (N 2 



128N r 



N c A g 
lN f A q 



'-(2N C + N f ) 



Ag - Ag 
A q - Ag 



N c A g 

l^fAq 



(6-30) 



In the absence of turbulent fields, this array of equations reduces to eqs. (6.8), (6.9) 
of Ref. 6). The result for Scenario B is obtained by means of the substitution (6-27). 

It is useful to introduce a concise vector notation for these equations. With the 
help of the two- vector A = (A g ,A q ), (6-30) can be written in the symbolic form 



(o A + a c )A = r , 



with 



r = 



3205)/ iV 2 -l 
3vr 2 V f N * N f 
32C(4) g 2 {B 2 )rT g 



5vr 2 



T 3 



N c 
7 



N f 



(6-31) 

(6-32a) 
(6-32b) 



20 



M. Asakawa, S. A. Bass, and B. Mutter 



7T(^c-l) 4, 



2,_ . N c 



9" ^' "\ ^ 
n 2 N f (N 2 - 1) / 1 -1 



128iY r \ -1 1 



(6-32c) 



The equation (6-29) for the shear viscosity then takes the form of a scalar product 
between the vectors r and A: 

»7= ^ r ' A= \ r ■ (aA + ac)" 1 (6-33) 

where we used (6-31) to express A in terms of r and the matrices a\ and ac- 

The inversely additive property of the contributions of turbulent color fields and 
parton collisions to the total shear viscosity is apparent in (6-33). Its origin lies in 
the additivity of the relaxation rates due to different processes and allows us to write 
the total shear viscosity as an inverse sum of anomalous and collisional viscosity: 

r?- 1 « Va 1 + rjrf , (6-34) 

with 

3 3 

VA = ^r-a A 1 -r, rjc = j r ■ a^ 1 ■ r . (6-35) 

Equation (6-34) implies that, whichever contribution to rj is smaller, dominates the 
overall shear viscosity The anomalous viscosity dominates when t]a < rjc- Because 
r/A grows with a smaller power of g than rjc, the anomalous viscosity dominates for 
sufficiently weak coupling, and we thus have rj ~ 77 a for g — > 0. 

6.6. Estimate of the anomalous viscosity 

Equation (6-21) shows that r/A decreases with increasing strength of the turbulent 
fields. Since the amplitude of these fields grows with the magnitude of the momentum 
anisotropy, a large anisotropy will result in a small value of r/A- The anomalous 
mechanism thus exhibits a stable equilibrium in which the shear viscosity regulates 
itself: The momentum anisotropy grows with r/A, but an increased anisotropy tends 
to reduce the anomalous viscosity. This leads to a self-consistency condition which 
determines r/A- 

In order to proceed further, we need to explore the dependence of the turbulent 
fields on the anisotropy of the momentum distribution of the partons. In particular, 
we must know how large (B 2 ) and r m are. The coherent color magnetic fields are 
only generated by the plasma instability when the momentum distribution of partons 
in the quark-gluon plasma is deformed due to the collective expansion. Analytical 
studies have shown that the instability always occurs when the momentum distri- 
bution is anisotropic. We also know from these studies how the growth rate of the 
instability depends on the anisotropy of the momentum distribution, but there are 
no published systematic studies that show how the "saturation" level of the coher- 
ent field energy depends on the anisotropy. We will therefore rely on some heuristic 
arguments for the needed dependences. 
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The study by Romatschke and Strickland 9 ) uses the following parametrization 
of the anisotropic momentum distribution: 



/(p) = fo (Vp 2 +t(vn) 2 ) « Hp) - ^^/o(l ± So) • (6-36) 



Choosing h = e z and subtracting the trace, this corresponds to a perturbation of 
the equilibrium distribution: 

«P> = -db( ri -T)- (6-37) 
Comparing with (5-12) this establishes the connection 

C = 2Z\^. (6-38) 

The relative anisotropy of the stress tensor is given by 

5T rT ST in , oT zz 8 . , . 

2 # = 2 # = = 15^ (6 ' 39) 

^11 J ?/J/ J 22 

Comparing with (5T4) we obtain a relation between ry and £, which takes the form 
(for a massless parton gas): 

5= i^v !i | =10 2jvu| (fj . 40) 

The central point of our argument is that the average collective color field energy 
is a function of the momentum anisotropy. For Scenario A (turbulent color-magnetic 
fields only) this means: (B 2 ) = £>(£). We do not know this function in detail, but 
we know that 6(0) = 0, because no instability exists in the absence of a momentum 
anisotropy The simplest ansatz is a power law, which we will write in the form 

g 2 (B 2 ) = b g 4 T 4 C (6-41) 

with an as yet unknown power n. Following the discussion before equation (3-7) we 
set n = 2.*) The memory time r m can be determined by one of two mechanisms. If 
the plasma particles move faster than the coherent fields propagate, r m will be set 
by the spatial coherence length of the coherent fields. This coherence length is given 
by the wave length of the maximally unstable mode, which is of the order of the 
Debye length:**) 

r m = d r 1/2 (gT)- 1 . (6-42) 



*•* This differs from the assumption made in Ref. 5), where we assumed n = 1. 
**' In the opposite situation, when the color fields evolve rapidly compared with the motion of the 
plasma particles, the memory time is determined by the decoherence time of the fields interacting 
with the particles and among themselves via the nonlinearities of the Yang-Mills equations. The 
case, in which the back reaction of the particle motion on the unstable field modes determines the 
decoherence time has been treated by Dupree 73 - 1 and by Abe and Niu. 77 '' 78 ' 
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The decoherence time and its dependence on £ can, in principle, be determined 
from simulations of the classical Yang-Mills equations. In the absence of such a 
determination, we shall assume that r m is given by (6-42). 

We are now ready to derive the self-consistency condition for the anomalous vis- 
cosity in the limit when the collisional viscosity can be neglected. We are interested 
in the late time, steady state situation where the expansion drives the momentum 
distribution into a slightly oblate anisotropy along the z-axis. As discussed before, 
this leads to the growth of unstable mean field modes, which determines the viscosity 
ija which, in turn, controls the size of the momentum space anisotropy. The feed- 
back loop is a stable one, because a large anisotropy leads to large saturation levels 
of the fields, which reduce the viscosity and thus limit the size of the momentum 
anisotropy. 

One technical complication of the limit of vanishing collisional viscosity is that 
the momentum distributions of quarks and gluons then attain different anisotropics, 
because color charges differ and thus their interaction with the turbulent fields is dif- 
ferent. This implies that the presently unknown constants 60 an d do receive different 
contributions from quarks and gluons, and that we should, in principle, distinguish 
between £ ff and £ q . In order to avoid unnecessay distractions from our argument, and 
because the numerical constants are not known anyway, we will only consider the 
gluon contribution to the anomalous viscosity in the following. The generalization 
to include quarks is straightforward and does not change the functional dependence 
of the result on g, T, and |V«|. 

Combining (6-41) and (6-42) into a single scaling relation and replacing £ with 
rj = t/a, we obtain: 



We now insert this result into the expressions (6 -22a) for the anomalous viscosity 
and obtain the self-consistency relation: 



where we have combined all numerical constants into a single factor cq . We note once 
again that this relation neglects contributions from the collisional shear viscosity rye 
and does not distinguish between the contributions from quarks and gluons. We 
also note that the power on the right-hand side differs slightly from that in Ref. 5), 
because of our different scaling assumptions in (6-41) and (6-42). Resolving (6-44) 
for tja/s finally yields the desired expression for the self-consistent anomalous shear 
viscosity: 



g 2 (B 2 ) r m = b Q d Q {gTfe 12 = b d (gTf ( 




(6-43) 




(6-44) 




(6-45) 



Several things are noteworthy about this result. First, if the memory time r m is 
longer than our estimate (6-42), the value of t/a decreases. Second, the dependence 
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Fig. 1. Schematic representation of the dependence of the collisional and anomalous viscous stress 
STik on the velocity gradient. The collisional viscous stress is shown by the linearly rising 
black line; the anomalous viscous stress is shown by the (colored) curved lines for three scaling 
exponents of the turbulent color field energy (n = 1.5, 2, 2.5). Solid lines indicate the dominant 
source of viscous stress in different regions of the scaled velocity gradient |Vu|/T. 



on the gauge coupling (~ g & ^) is parametrically much weaker than that of the 
collisional viscosity (~ (g 4 ln(l/g)) -1 ). Thus, for early times r and weak coupling 
g <C 1, the anomalous viscosity will be much smaller than the collisional viscosity 
and thus dominate the total viscosity according to (6-34). 

It is instructive to consider the role of the anomalous viscosity in the hydrody- 
namical equations. As (5-14) shows, the viscous contribution to the stress tensor is 
proportional to ry|Vtt|. For the collisional shear viscosity which is, in first approxi- 
mation, independent of the magnitude of the velocity gradient; this implies that ST,-^ 
grows linearly with |Vu|. The anomalous shear viscosity (6-45), on the other hand, 
is a decreasing function of the velocity gradient; its contribution to the stress tensor 
grows like |Vu| 2 / 5 for our scaling assumptions. The unusual dependence of t/a on 
\Vu\ certainly justifies the term "anomalous viscosity". 

The different dependence of the collisional and the anomalous viscous stress 
tensor on the velocity gradient is shown schematically in Fig. 1. For very small 
gradients the linear dependence of the collisional viscous stress tensor dominates, 
but for larger velocity gradients the lower power associated with the anomalous 
shear viscosity will assert its dominance. The precise location of the cross-over 
between the two domains depends on the value of the numerical constant cq, but we 
can deduce from (6-45) that the cross-over point shifts to lower values of |Vu| with 
decreasing coupling constant g. We also show in the figure the effect of choosing a 
different power n in the scaling law (6-41) for the energy density of the turbulent 
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color fields. Finally, we note that the decreasing dependence of t/a on |Vu| implies 
that the shear viscosity can, in principle, fall below the Kovtun-Son-Starinets (KSS) 
bound rj/s = (47r) _1 . 

§7. Conclusions and Outlook 

In this paper we have presented details of our derivation of the anomalous vis- 
cosity in an anisotropically expanding quark-gluon plasma, which arises from in- 
teractions of thermal partons with dynamically generated color fields. In the weak 
coupling limit or for large velocity gradients, the anomalous viscosity is much smaller 
than the viscosity due to collisions among thermal partons. By reducing the shear 
viscosity of a weakly coupled, but expanding quark-gluon plasma, this mechanism 
could possibly explain the observations of the RHIC experiments without the as- 
sumption of a strongly coupled plasma state. 

Due to the self-consistency condition, the anomalous shear viscosity itself is in- 
versely dependent on the expansion rate of the plasma. This means that the viscous 
term in the hydrodynamic equation does not depend on the velocity gradient lin- 
early, but sublinearly. This unusual dependence on the velocity gradient justifies 
the term "anomalous" viscosity. It also implies that the response of the medium to 
the expansion is nonlinear, in fact less than linear. The usual method to relate the 
transport coefficient to a correlation function (Kubo formula) that can be measured 
in full equilibrium is therefore not applicable. Even if the appropriate real-time cor- 
relation function could be determined from Euclidean lattice QCD simulations, it 
would fail to describe the anomalous viscosity. More sophisticated methods, includ- 
ing a treatment of the color instabilities in the expanding quark-gluon plasma, will 
be necessary. 

The presence of anomalous viscosity for a rapidly expanding anisotropic quark- 
gluon plasma has important consequences for the early Universe. The expansion 
of the quark-gluon plasma in the early Universe is relatively slow and thermaliza- 
tion is therefore always maintained. Thus, the quark-gluon plasma in the early 
Universe does not generally possess an anomalous shear viscosity. However, when 
off-equilibrium processes induce a strong local anisotropy, the anomalous viscosity 
can become important in determining the total viscosity. This could have been the 
case during the reheating period after inflation and the electroweak phase transition, 
and may have affected the production of fluctuations in the Universe, baryogenesis, 
etc. 

The approach described here can be applied to other transport properties of an 
expanding, turbulent quark-gluon plasma. Examples are the coefficient q of radiative 
energy loss of an energetic parton, 95 ) the absorption coefficient for bound states of 
heavy quarks (e. g. 3/ip), and flavor equilibration rates. For the rate of strange quark 
pair production, we note that a quark pair can be produced in the presence of a mean 
color field by a single gluon. This process is well known in quantum electrodynamics, 
where a single photon can produce an electron-positron pair in the presence of a 
strong electromagnetic field, such as the Coulomb field of a nucleus. In QED the 
pair production rate grows like Z 2 , where Ze is the nuclear charge; in the case of 
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the quark-gluon plasma one would expect the rate to grow like g 2 (B 2 } and hence be 
proportional to the expansion rate. A rapid expansion would, therefore, accelerate 
the rate of light quark pair production. Such an effect may, indeed, contribute 
to the very high rate of quark pair production seen in numerical solutions of the 
Dirac equation in longitudinally expanding gluon fields in the color glass condensate 
picture. 96 ) 

Turbulent color fields can also contribute to the dissociation of heavy quark- 
antiquark bound states, because the two constituents will be deflected in oppo- 
site directions. This effect has been considered for randomly oriented color-electric 
fields, 97 ^ but not for color-magnetic fields, which should have a similar effect, albeit 
suppressed by a factor (v/c) 2 , where v is the velocity of the charmonium. Turbulent 
color fields will also contribute to the diffusion coefficient for heavy quarks in the 
plasma. Finally, turbulent color fields can influence the trajectories of the partons 
contained in a jet created by a hard scattered quark or gluon. Because the fields cre- 
ated by the color instabilities in an expanding medium are polarized transversely to 
the expansion direction, the color magnetic fields will preferentially deflect the out- 
going partons in the longitudinal direction and thus cause a longitudinal broadening 
of the jet cone. Such an effect has been been observed at RHIC. 98 ) 

In a forthcoming publication we shall further investigate the phenomenological 
consequences of the anomalous viscosity: in particular we shall focus on jet energy- 
loss in turbulent color fields and address the question whether the vast array of 
current data at RHIC can be understood in the framework of a weakly coupled QGP 
with anomalous transport properties. 
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Appendix A 

Alternative Derivation of the Vlasov Diffusion Term 



In this section, we present a different derivation of the Vlasov diffusion term 
(4-14), based on the evolution of the extended phase-space distribution f(r,p,Q,t), 
which is given by 81 ) 



^ + 9f abc Q a A b ^ + 9 Q a (£° + v x • V p 



f(r,p,Q,t)=0. (A- 1) 



Here we have suppressed any unnecessary vector notation and omitted the colli- 
sion term. Our alternative derivation of the diffusion term starts from the integral 
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representation for /: 

f(r,p,Q,t) = j dr dpodQo8{r-f{t))8(p-p(t))S(Q-Q(t))f(ro,po,Q ,0), (A-2) 
with 

df _ p 

l~t =v = E- p > (A ' 3a) 

f| = g Q a (S a {r) + vx B a {f)) , (A-3b) 
^§ = gf abc Q h A c X- (A-3c) 

and ro = f(0), p(0) = po, Q(0) = Qq. Here f(t) denotes the trajectory of a plasma 
particle, which is found at position ro with momentum po and color charge Qq at 
time t = 0. The evolving phase space distribution is given by the sum over all the 
trajectories of these (test) particles. When we act on (A-2) with the Vlasov drift 
operator (d t + v • V r ) and use the definition (A-3) of v, we obtain: 

{jt + v ■ Vr + afabcQaA y^) ^ q > *) 

= - j dr dp 5(r - r) -VpSip - p) f(r o ,p o ,Qo,0), (A-4) 

where we dropped the explicit notation of the time-dependence of f, p, and Q. Since 
p is time independent for vanishing color field, it makes sense to expand the delta 
function 5{p — p) around (p — po) for weak fields (or short times). It is important 
here to remember that (strong) plasma turbulence does not imply strong fields, just 
a spectral field distribution without phase correlations on all inverse length scales k. 
The expansion is: 

5( P -p)= 5(p - po) - Ap(t) ■ V p 5(p - po) + • • • , (A-5) 

where Ap(t) = p(t) — po, implying d{Ap)/dt = dp/dt. The right-hand side of (A-4) 
then becomes: 

- J dr dpo 5(r - f) . V p [1 - Ap(t) • V p ] S(p - p )f(r ,p , Qo,0) +• ■ ■ . (A-6) 

The expression now contains a term quadratic in Ap, and it makes sense to take an 
ensemble average over the color fields. We also integrate by parts with respect to 
V p and obtain: 

( » + v . v , + *) ~ < v > • Avit)i -1T ■ v ' ') 

= Vj-D-Vpi/). (A-7) 
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In the last step we have introduced the diffusion coefficient D and assumed that we 
can factorize the ensemble average. Explicitly, the diffusion coefficient is given by 

(i,3 = x,y,z): 

= (Apiit)^) = g? ^dt'(Q a (t')F^f(t'),t')Q\t)F^f(t),t)), (A-8) 

where 

F?(r(t),t) = £i(.m,t) + (Ht) x B a (f (t), t)\ . (A-9) 

In the local rest frame of the medium the electric and magnetic components of the 
color field can be expressed as (using the convention e 0123 = 1): 

£? = F? v u v ; B? = e iX ^F aX ^u v . (A- 10) 

If we now argue that the correlation time/length for the color fields is short in 
comparison with the temporal change of the velocity of a plasma particle, we can 
take v(t) = v(t') = v out of the average and are left with the autocorrelation function 
of the color fields and color charges along a typical particle trajectory. 

The time evolution of the color charge Q a (t) is given by the solution of the last 
equation (A- 3): 

Q a (t') = P exp (j^ ) f abc A b ^dx^j Q c (t) = U ac (f(t'),f(t))Q c (t) . (A-ll) 
Inserting this solution into (A-8), we obtain 

D '' ^ (4^) 

= g 2 £ alt' (Q%t)Q\t)F?{f{t%t')UUf{t% f(t))I*(f(t), t)) . (A-12) 

If we now assume that the distribution of the color charges of partons at a given 
time t is random and independent of the color fields, (Q a Q b ) = (iV 2 — l)" 1 ^^, we 
recover the expression (4-14) for the Vlasov diffusion term with the field correlators 
(4-11). 

Finally, in order to make contact with the diffusive transport equation (5-1), we 
need to integrate (A- 7) over the color charges Q to obtain an equation for 

f(r,p,t) = J dQ(f(r,p,Q,t)). (A-13) 

Partial integration with respect to Q shows that the third term on the left-hand side 
of (A-7) vanishes in view of the relation dQ a /dQ c = 5 ac and the antisymmetry of 

fabc' 

We end this section with a remark concerning abelian plasmas. Our derivation 
of the diffusive Vlasov term was motivated by the insight that the color instabilities 
of an anisotropic nonabelian plasma saturate under the action of the nonlinearities 
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of the Yang-Mills equation.")' 10 °) This mechanism is, of course, absent in abelian 
plasmas. The saturation of the instablilities is then caused by the back reaction of 
the growing soft field modes on the particle distribution. 77 ) In order to address 
this situation, we follow Dupree's argument 73 ) that the ensemble average in (A-8) 
should be taken over Fourier components of the color field, because these are the 
slowly varying variables. This is standard practice when dealing with an ensemble 
of waves. It certainly makes sense for electromagnetic plasmas, where nonrelativistic 
particles move in fields that propagate at the speed of light. In the case of the 
quark-gluon plasma, however, the situation is reversed: the thermal partons move 
with the speed of light, but the soft color fields obey a dispersion relation with a 
slower propagation speed. Thus it is doubtful whether a similar reasoning would 
make sense. Nevertheless, for the sake of interest, we outline Dupree's approach (for 
color-magnetic fields only). One writes 

B a (f(t),t) = Y J B a {k)e- i ^ t+ik < t \ (A- 14) 

k 

where B(k) are the Fourier components of the field. One can then pull the slow 
variables B(k) out of the time integral in (A-8). After factorizing the ensemble 
average one obtains 

V p ■ D ■ V p = -<?<?<?Y,W( k )&A- k )) L ® L T (J* dlfe^-V^W^®^ . 

(A-15) 

The time integral can be interpreted as an autocorrelation or memory time for the 
action of the magnetic fields on the particles. In Dupree's approach, the value of 
the time integral is governed by the effect of the turbulent magnetic fields on the 
particle trajectory, and can be shown to satisfy a self-consistency condition. We will 
not pursue this approach further here and refer the interested reader to Abe and 
Niu's work. 78 ) 

Appendix B 

Color Instabilities near Equilibrium 

For convenience, we here state the results for the growth rate of unstable quark- 
gluon plasma modes in a background distribution of quasi-thermal partons whose 
momentum distribution is only slightly perturbed away from equilibrium. We follow 
the notation and derivation of Romatschke and Strickland. 9 ) 

We assume that the momentum distribution in the rest frame of the medium 
can be written as 

/(p) = fo (Vp 2 + £( P -n) 2 ) , (B-l) 

where n is a unit vector defining the orientation of the anisotropy and £ <C 1. Wc 
denote the angle between the wave vector k of the considered field mode and n 
by 9: cos 9 = k • n, where k is the unit vector in the direction of k. For a given 
temperature T, the polarization function of the gluon field is then a function of 
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the variables u>, k = |k|, and 9. Because the parton distribution violates spherical 
symmetry, there are four different components of the polarization tensor, which can 
be expressed in terms of the functions a, f3, 7, 5. The gluon propagator in medium 
can be decomposed with the help of the projectors i, j = 1, 2, 3, 



A - — (S- ■ 



Bj 



kjkj , 



k%kj , 



Di 



(B-2a) 
(B-2b) 
(B-2c) 
(B-2d) 



/% 3 — kjjij -\- fiikj , 

where h = n — k(k • n). The resulting expression for the gluon propagator is 

Aij(k, u, 6) = A K {Aij - dj) + A G [{k 2 -Lo 2 + a + 7)^ + (J3- uj 2 )C %1 - 5D tj ] (B-3) 
with 

A^[k,w,0) = k 2 - to 2 + a, 



&£{k,u,0) = (k z - u z + a + - u z ) -ri z 8 z . 



(B-4a) 
(B-4b) 



For small values of £ the functions a, (3, 7, 5 are given by 



- iW) y(5 cos 2 ^ - 1) - ^ cos 2 # 



/? = z 2 77 L (z)+£z 2 



+i/7 T (, z )(3cos 2 0-l) 



^-z 2 i7 L (z))(3cos 2 ^ L) 



7 
<5 



+i7 L (z) (2 cos 2 6» - 

e(^ T (z)-^) {z 2 -i)^ 2 e, 

e( 4 z 2 ^-77 T (z)(l-4z 2 )) cos#, 



1) 



where z = u/k and 



n T (z 

n L (z) = mjj 



D J2 



2 

,2 



z 



, z + 1 
. m 

2zV z-l 



£ini±i-i 

2 z-l 



(B-5a) 

(B-5b) 
(B-5c) 

(B-5d) 

(B-6a) 
(B-6b) 



are the usual expressions for the transverse and longitudinal gluon polarization func- 
tions in thermal equilibrium. 
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